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Abstract 

We present a study of patterns, formed in drying drops of aqueous gelatin solution 
containing sodium sulphate. The patterns are highly complex, consisting of a hierar¬ 
chical sequence of rings which form concentric bands as well as dendritic crystalline 
aggregates. When the preparation of the complex fluid is done by mixing an aqueous 
solution of the salt with an aqueous solution of gelatin prepared separately, another 
feature is observed in the pattern on the dried out drop. This is a viscous hngering 
pattern, superposed on the series of rings. We try to explain the origin of these two 
features from a simple physical approach. 
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1 Introduction 


The study of crystallization patterns of salts in drying droplets shows a wide va¬ 
riety of intricate and beautiful patterns^’^. Sodium sulphate and sodium chloride, 
in high viscosity solutions are popular candidates for such studies. Crystal forma¬ 
tion on macroscopic scales to microscopic dendritic growth have been reported^, 
ring formations like the well-known coffee stain‘d and multiple rings have also been 
observed^. Many different aspects of drying droplets are currently under study®’^. 

Sodium sulphate is a constituent of Portland cement and usually crystallizes as 
thenardite (Na 2 S 04 ) or mirabilite (Na 2 SO 4 . 10 H 2 O)® depending on ambient con¬ 
ditions and water content. Several different polymorphic forms have been ob¬ 
served. Rodriguez-Navarro et al.^ further show that crystallization pressure devel¬ 
oped within porous materials when thenardite is formed, can cause more damage to 
building material and rocks than formation of mirabilite. A metastable form, sodium 
hepta-hydrate (Na 2 S 04 . 7 H 20 ) is also believed to be signihcant in this context®’®. So 
such studies are important from a practical point of view. 

If the salt solution in pure water is dried, a coffee-ring pattern‘d is observed. On 
highly hydrophillic surfaces, the salt forms dendritic patterns outside the drop ini¬ 
tially deposited®. This is due to a precursor wetting him surrounding the drop. The 
anhydrous salt crystallizes in rhombohedral form. If the salt is dissolved in a com¬ 
plex viscous huid e.g. containing a polymer, the dried drop shows more complicated 
features®. 

In the present paper we focus primarily on two interesting observations - (i) 


a hierarchical series of concentric rings whose distribution is similar to a devil’s 
staircase^^ and (ii) a viscous hngering pattern^^ superposed on the rings, which is 
observable under certain conditions. 


2 Materials and Methods 

The complex solution for droplet deposition was prepared by two slightly different 
methods and the results showed a noticeable difference in the pattern left on the 
substrate after the droplet dries. 

2.1 Method I 

0.142 g sodium sulphate powder was mixed with 50 ml. HPLC grade water, to get a 
solution of 0.02M concentration. 0.5 g gelatin powder was mixed with this solution 
and stirred at temperature 60°C until completely dissolved. 

2.2 Method II 

0.2 g of gelatin was added to 20 ml of HPLC grade water, this gives a concentration 
of 1% by weight. It was stirred in a magnetic stirrer at temperature of 60°C until 
it completely dissolves. 

Sodium sulphate solution in water was prepared by mixing a; g of salt in 20 
ml. water, where x = 0.028, 0.056, 0.085, 0.113 and 0.142 respectively. These two 
solutions were mixed in 1:1 ratio by volume. The molar concentration of the salt 


was varied from 0.01 to 0.05 M. 


A droplet of each solution of volume ~ 50/il was deposited using a micro-pipette 
on a clean glass surface and allowed to dry under ambient conditions. 

3 Results 

The appearance of the drops on drying and some relevant physical properties were 
measured. Results are given below. 

3.1 Physical properties: Surface tension, Viscosity and an¬ 
gle of contact measurement 

Surface tensions of (i) sodium sulphate in water, (ii) gelatin in water and (iii) the 
fluid containing both salt and gelatin, prepared by method II, were measured by 
KRUSS Tensiometer K9 at Physics Department, Jadavpur University. Viscosity 
of the above three solutions was measured by Brookheld DV-III Ultra viscometer 
at Gurudas College, Kolkata. The angles of contact of sample (iii) on glass was 
measured to be approximately 18 ±2°. All results are given in Table ([T]) 

3.2 Imaging 

Gross features of the pattern on the dried drop were visible to the naked eye, but 
new features were revealed as the resolution of observation changed. 


3.2.1 Low resolution imaging 


Droplets dried on glass show a pattern which depends on the method of preparation 
of the sample. The dried droplets were photographed using Nikon CoolPix L120. Fig 
([I]) shows typical patterns generated. When the aqueous gelatin solution is prepared 
and dry powdered salt mixed with it, the patterns are as shown in Fig (II])(a) for 
salt concentration 0.01 M. A series of concentric rings are formed, surrounded by 
a transparent annular region. When the salt concentration is increased to 0.03 M, 
the pattern looks more or less similar but the transparent border is much reduced 
in width (hg[T] (b)). 

However, when the sample is prepared with the same constituents, but following 
a slightly different procedure, a new feature is observed in the pattern. Now the 
gelatin-water and sodium sulphate-water solutions are prepared separately and then 
mixed. In this case a radial hngering pattern appears, which is superposed on the 
rings (hg ([I])(c) and (d)). The hngering pattern is seen for low salt concentrations 
~ 0.02M and it disappears when salt concentration is increased to ~ 0.04M or 
higher. 

3.2.2 Medium resolution imaging 

The rings shown in hg(IT]) are actually bands of several closely spaced thinner rings. 
The rings are not equispaced, but the spacing has a certain order. This hue structure 
is seen at higher resolution. Complex dendritic structures growing from the rings 
are also visible, when the drops are viewed under a microscope (Leica DM 750). 
Fig ([2]) shows the details of the pattern, for different regions of the droplet. What 


appeared to be a single rings in fig ([T]) is seen to be sets of rings with dendritic 
crystals growing in between them. The fine structure of the different regions is 
shown in Fig ([2]) where relative positions of the different micro-graphs are explained 
with respect to a schematic diagram of the whole pattern. 

3.2.3 Scanning Electron Microscopy (SEM) 

FESEM imaging has been done using SEM facility, configuration no. QUO-35357- 
0614 (funded by FIST-2,DST Government of India, at the Physics Department, 
Jadavpur University) to see the structure of the dendrites down to nano scales. The 
structures for increasing resolution are shown in fig ([3]). Diamond shaped crystallites 
of sodium sulphate can be discerned at highest resolution as shown in fig ([3](d)). 

4 Analysis of patterns 

The ring patterns have a typical hierarchical structure which looks self-similar for 
a certain range of length-scales. To analyse the pattern, the following procedure is 
adopted. The radial position of each ring with reference to the centre of the pattern 
is tabulated. Let us number the rings consecutively from the central reference point 
and denote the radius of the nth ring by i?„. We plot n on the Y-axis against on 
the X-axis, as shown in fig (jl]). If the points are joined to the Y-axis by horizontal 
lines, such that all points where Rn < Rn+i are on the same ‘step’ corresponding to 
Y = n, we get a staircase-like structure. It appears that as we observe the pattern 
at higher and higher resolutions new rings appear, that is new steps appear between 


the existing ones. So this pattern resembles a Devil’s staircase where ideally an 
inhnite nnmber of steps exist between any two steps. For real systems of conrse, 
there is a hnite cut-off. 

It is seen from £g (jl]) that there are a few wide steps interspersed by many shorter 
steps. We plot a cumulative distribution of step lengths by counting the number 
Ni of steps with length > I and plotting Ni against 1. The resulting curve can be 
approximately fit by a logarithmic function of the form 

Y = -A- BlnX (1) 

as shown in £g ([S]). 

4.1 Hierarchical ring pattern: a simple model 

The hierarchical pattern of rings is a notable feature of the pattern and we offer a 
simple explanation for the process. Kaya et al.^ have calculated the observed regular 
ring spacings for the system sodium poly(styrene sulfonate) with sodium chloride. 
Their approach is as follows: the ‘coffee ring’ at the outer periphery of the droplet 
is formed due to the highest evaporation rate at the outer boundary of the drop, 
which draws fluid carrying the solute towards the boundary. The solute particles get 
pinned at the boundary to form the so called ‘coffee-ring’. Kaya et ah suggest that 
the process may not stop here. The high concentration of solute at the boundary 
drives a diffusion flux inward, this superposes the outward advective flux and gives 
rise to non-monotonic concentration contour in the region nearer the droplet centre 


from the coffee ring. Kaya et al.^ calculate the shape of this concentration contour 
and hnd a maximum, which exceeds the critical solute concentration for deposition. 
This leads to a second ring of smaller radius than the outer one and concentric with 
it. This process repeats to form a third ring, fourth ring and so on. In the model 
of Kaya et ah only the largest ring formed, participates in the diffusion process and 
the resulting rings are approximately equally spaced. 

We offer a simple heuristic extension of the above scenario. We assume that 
when the fluid is complex and highly viscous, like gelatin, the outer rings do not 
dry out immediately and so they continue to participate in diffusion. We suppose 
that all, or at least several rings are active simultaneously. We show by a simplihed 
argument that this leads to a system of rings whose spacings closely resemble the 
real pattern generated in our drying experiments. 

From the concentration contour calculated by Kaya et al.®, we see that the 
contour is asymmetric with a higher gradient outward (towards the ring boundary) 
and a lower gradient inward. So we may assume faster diffusion outward and slower 
inward and as a result when the two concentration fluxes superpose, the maximum 
of the resultant contour will be shifted towards the larger ring, instead of forming 
exactly in the middle. Let us start with an outer ring at radius ri and a subsequent 
one at r 2 (see hg[6]), where the radii are related by a constant ratio b = r^/ri. We 
assume the next ring to form at ra, due to the superposing diffusion fluxes from 
rings 1 and 2, maintaining the same ratio. So 



this process follows for successive rings which form closer and closer to ri, for 


the nth and (n + l)th such rings we have then 


rn+l - Tn 

ri - Tn 


(3) 


This allows us to write a recurrence relation for the radius of successive rings as 


Tn+i = rn{l-b) + hri (4) 

for n > 1. 

We can thus generate a sequence from two initial parameters b and ri Similar 
such sequences can be written taking the initial ring radius as r 2 , with n > 2 and 
so on with the ratio b replaced by other constants say c, d etc. 

Then the next sequence is given by 

rm+i = rm{^ - c) + cr 2 (5) 

for m > 2. 

Plotting the positions of several such series we get a ring pattern like £g (j7j). 
where b is varied successively from 0.6 to 0.9 (which is close to the experimental 
results). It can be seen that the resulting pattern of rings qualitatively reproduces 
the experimental patterns (hgHj). In the real system, dendrites grow in the space 
between the adjacent rings, as shown with increasing magnihcation in SEM images 
(fig ED . We calculate the cumulative number of steps Ni greater than a length / and 


plot it against /, as done for the experimental ring sequence. The results can be fit 



approximately to a logarithmic plot as shown in fig (|8]). So this very simple pictnre 
gives a qnalitative description of the experimental patterns. 

4.2 Viscous fingering pattern 

The low resolution images shown in fig([I]) show fingering patterns superposed on 
the rings, when the solution is prepared by Method II. It is well known that viscous 
fingering patterns are formed when a low viscosity fluid, at a higher pressure is 
forced into a higher viscosity fluid, kept at a lower pressure^^. The two fluids should 
be immiscible. 

The formation of the fingering pattern in the drying droplet may have the fol¬ 
lowing origin. When Method II is employed the fluid (A) - a solution of gelatin in 
water is mixed with the solution (B) which is a solution of the salt in water. As the 
droplet dries, there is a tendency for the solutions to separate, this is seen by the 
formation of the transparent borders in fig([I]), which are mostly gelatin. Viscosity 
measurement shows that, before drying starts the salt solution has a lower viscosity 
(~ Icp) compared to the gelatin solution, which has viscosity (~ 1.6cp). As drying 
starts, a pressure gradient develops with a lower pressure at the drop periphery, due 
to higher evaporation rate here. So the condition for viscous fingering is satisfied 
and patterns form. Surface tension measurements show that the two fluids have 
different surface tensions and are therefore immiscible. Though the fluids appear 
to form a homogeneous solution when mixed together, on a microscopic level they 
retain their identity as in an emulsion. 

At higher salt concentrations, salt becomes more uniformly distributed, there is 


Sample 

Surface Tension 

Viscosity 


(mN/m) 

(cp) 

Gelatin solution 

39.0 

2.06 

Na 2 S 04 solution 

70.0 

1.65 

Mixture of Gelatin and Na 2 S 04 solution 

40.6 

2.15 


Table 1: Measured physical properties of the experimental solutions at 30°C and 
relative humidity 52% 

no gelatin border and the hngering pattern disappears. Following the preparation 
Method I also shows no hngering. In this case two different huids are not mixed, as 
the gelatin is added in powder form to the salt solution, so viscous hngering does 
not occur. 


5 Discussion 

The study of drying drops has become a highly popular topic in recent times^^’^’®’^^. 
The beauty of the fascinating patterns formed is not the only reason for this widespread 
interest. Medical diagnostics^^ promises to be an area where patterns of dried bi¬ 
ological huids will be of practical use. Other practical applications are in ink-jet 
printing and drying of thin coatings. 

We discuss here two aspects of these patterns which to our knowledge, have 
not yet been reported or explained. Though concentric rings have been reported in 
drying droplets, the typical pattern we observe has not yet attracted much attention. 






The simultaneous presence of viscous fingering and ring-dendrite patterns appears 
to be a new observation as well. The models we propose here are very elementary 
and need to be developed for a complete understanding of these phenomena. We 
hope to pursue these studies and report work on more detailed models in future. 
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Figure 1: The appearance of the dried droplets are shown, the graduations on the 
scales below are in mm. (a) and (b) are prepared by Method-I, while (c) and (d) by 
Method-II show viscons hngering. (a) is for salt concentration 0.01 M and (b) for 
0.03 M. (c) is for 0.02 M solntion of 50 /r/ volnme, while (d) is for 0.02 M solntion 


of 40 fil volnme. 




Figure 2: Images of different regions of the pattern under microscope are shown. 
The scale bars in the hgure panels are all 100 pm. Hierarchical ring patterns with 
dendritic crystal growth in between are observed. 



Figure 3: SEM images taken at successively increasing resolution from (a) to (d) are 
shown. High resolution images of the dendrites show sodium sulphate crystallites. 





Figure 4: The measured radii of the successively increasing rings are shown as a 
‘staircase’ pattern. The inset shows a zoomed up portion, demonstrating the self¬ 
similarity. 



Figure 5: The number of rings Ni with length > I are plotted against 1. A logarithmic 


ht is shown in red. 











Figure 6: Schematic diagram showing successive ring formation. The ratio r 2 /ri = b. 
The contours below are approximate outlines of the variation in salt concentration 
with distance from the center (after Kaya et al.^) and the arrows represent the 
opposing diffusion fluxes. A ring forms when the maximum of the contour exceeds 


the critical salt concentrstion. 









Figure 7: Radii of rings calculated from eqs (jl]and[5]) are shown, here the ratio b is 
successively varied between 0.6 to 0.9. 



Figure 8: The number of rings Ni with length > I from hgd?]) are plotted against 1. 


A logarithmic £t is shown in red. 











